Study of shape transitions in N ~ 90 isotopes 
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Abstract. We study the spherical to prolate-deformed shape transition in 144 ~ 158 Sm and 146 - 160 Gc[ 
isotopes with modern calculations beyond the mean field with the Gogny D1S force. We compare 
the results with the shape-phase transition predicted by the collective hamiltonian model and with 
the experimental data. Our calculations do not support the existence of a first order phase transition 
in these isotopic chains in the viewpoint of the Bohr hamiltonian neither the interpretation of the 
nuclei N = 90 as critical points. 
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Recently, the study of the shape transitions in isotopic chains and its connection with 
possible quantum phase transitions has been an intensive field of research. The renewal 
interest about this topic is partly due to the development of analytic solutions for the 
critical points that correspond to phase transitions between spherical-prolate deformed 
(X(5)) and spherical- 7- soft (E(5)) systems in the context of the Bohr hamiltonian 03 
0. In particular, transitions between nuclei with vibrational spectra to nuclei with a 
rotational character have been observed experimentally in the N ~ 90 stable isotopes. 
Furthermore, the transitional nuclei 150 Nd, 152 Sm and 154 Gd have been proposed as 
the empirical realization of the critical points of first-order phase transitions due to the 
good agreement between the experimental data and the predictions given by the X(5) 
model (31 HI O. In this contribution we study the spherical-prolate deformed shape 
transition in the isotopic chains corresponding to 62Sm and 64Gd with state-of-the-art 
beyond mean field calculations using the Gogny D1S interaction. This study extends the 
analysis performed for the Nd isotopes in the Ref. [7] to the Sm and Gd nuclei. We will 
compare the theoretical results with the experimental data and also we will analyze the 
possible existence of phase transitions in this region of the nuclear chart. 
The description of these shape transitions as shape-phase transitions has been performed 
in the context of the Bohr collective hamiltonian. In particular, the spherical-prolate 
deformed transition is obtained from the following potential in the (/3 , 7) plane, deduced 
from the classical limit of the IBM model Q: 
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where N is the number of valence pairs -we take N = 10 as in Ref. [2 J- and £ (0 < £ < 1) 
is the continuous control parameter that is used to simulate the shape transition. In Figs, 
[jja-f) the evolution of this potential in the (j8, 7) plane as a function of the parameter £ 
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FIGURE 1. (a)-(f) Potential energy V(/3,y) (Eq. [I]) given by the collective model to describe the 
spherical to prolate deformed transition for different values of the control parameter £. Contour lines 
are plotted each 0.05 arbitrary units, the energy is set to zero at the minimum of each surface, (g) Value 
of the j3 deformation that corresponds to the minimum of (Eq.[T]i as a function of £. The inset shows the 
transitional region in detail. 



is shown. For £ = a spherical minimum is obtained and the potential is very stiff in the 
/3 direction. Increasing the value of £ the potential softens in the trajectory (/3, 7= 0°), 
emerging for £ = 0.025 a deformed minimum that coexists with the spherical one in a 
small range of £ values (Fig. [jjg)). For higher values of the control parameter a single 
deformed minimum is obtained and the /3 = point turns into a saddle point. Precisely 
the value of £ for which both minima are degenerated defines the critical point of the 
phase transition (see Fig. |4|a)). 

We now show the results of the beyond mean field calculations performed with the 
Gogny D1S force [6J. The microscopic equivalent to the Eq.[T]are the Potential Energy 
Surfaces (PES) |[T0l which are calculated for each nucleus evaluating the following 
expectation value: 

F N,z ( n v] = (^ z (/3,y)|^|^(/3,y)) 
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where \<t> N ' z (j3 , 7)) are HFB-type wave functions that are projected -before the variation- 
onto good particle numbers (N, Z) and constrained to the deformation parameters (/3 , 7) 
- see Refs. B71[8]| and references therein for more details. 

In Figs. [2j|3] the PES with number of neutrons between N = 82 — 96 are represented for 
Sm and Gd nuclei respectively. Similar results for Nd isotopes are obtained in Ref. 
and are not shown here. The three elements exhibit a similar general behavior. We ob- 
serve that increasing the number of neutrons a shape transition from spherical shapes 
-obtained for the semi-magic N = S2 nuclei- to well prolate deformed nuclei is given. 
However, the manner in which these transitions are produced is different to the results 
of the collective hamiltonian shown in Fig. [T] On the one hand, the microscopic calcu- 
lations reveal the importance of the triaxial degree of freedom in the transitional nuclei 



7=60° y= 60 ° Y= 60 ° Y= 60 ° 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.1 0.2 0.3 0.4 0.5 0.6 0.7 

P P P P 



FIGURE 2. (a)-(h) Particle Number Projected Potential Energy Surfaces E N ' Z (j3 , y) (Eq. [2} calculated 
with the Gogny D1S force for g^ 4 ~ 158 Sm isotopes. Contour lines represent a step of lMeV and the energy 
is set to zero at the minimum of each surface. 
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FIGURE 3. Same as Fig. gbut for ^ 6_160 Gd isotopes. 



N = 86, 88, 90 (Figs. [2]j3jc-e)). Here the PES are soft in the 7 direction in contrast to the 
description given by the collective hamiltonian, where the degeneracy is produced along 
the 7=0 trajectory. On the other hand, none of the PES obtained with the Gogny force 
shows the structure of a critical point in the sense that we do not obtain any flat potential 
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FIGURE 4. (a) Potential energy V Q3 , y = 0, 60) (Eq. [T]l in the critical point as a function of j3 . Vertical 
lines mark the position of the two minima, (b)-(d) Particle number projected potential energy surfaces 
along the axial trajectory calculated with the Gogny D1S force for 144 ~ 156 Sm and 146 ~ 158 Gd isotopes 
respectively. The zero of the energy is chosen at j3 = 0. 



in the (/3, 7 = 0) direction with two degenerated coexisting minima, one of them spheri- 
cal and the other one prolate deformed. This characteristic can be clearly seen in Fig.|4} 
where the PES along the axial direction are represented for Sm (Fig.[4jb)) and Gd (Fig. 
4 c)) nuclei as well as the collective potential that corresponds to the critical point in Eq. 
T (Fig.|4fa)). We observe that spherical minima are obtained only for the shell closure 
N = S2 nuclei. For the rest of isotopes we obtain oblate and prolate deformed minima 
where both the value of the /3 for which those minima appear and also the height of 
the barrier between them are increasing with higher values of the number of neutrons. 
Obviously, most of the oblate minima are indeed saddle points in the (/3,y) plane as it 
has been shown before. 

In order to obtain theoretical predictions not only for the PES but also for the spec- 
troscopy we carry out calculations within the framework of the Generator Coordi- 
nate Method (GCM) with the wave functions ^^(jS, y = 0,60)) also projected onto 
good angular momentum /. Hence, the resulting wave function can be expressed as: 
|vpiV,z,y,cT^ = j f^J,G(p)pJ\<$>N,z^ d p where pj is the angular momentum projec- 
tion operator and the coefficients f N > z ' J - a and the energies of the collective states 
E N ' Z ' J -° are found by solving the Hill-Wheeler-Griffin equation ||T0l 

J (^ N < z \P,P')-E N ^ J > ^ N ^\P,p ! ))f N > z > J >°(p ! W (3) 

being J^ NZ,J (P,P') and ,sV N,z,J (P,P') the projected hamiltonian and norm overlaps 
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FIGURE 5. Ratio between the energy of the first 4 + and 2 + excited states for (a) Sm and (b) Gd 
isotopes. The dash line is the X(5) prediction. 



respectively IfTOl . The calculations are performed with intrinsic axially symmetric wave 
functions due to the huge computational cost of the full triaxial case. The first micro- 
scopic calculations performed to study the shape transitions of the 148 - 15 °- 152 Nd iso- 
topes using a method very similar to the described above were carried out by Niksic and 
collaborators with a relativistic interaction [|9]|. 

The results obtained for the ground state band (c = 1) are compared with the experimen- 
tal data in Fig. [5} Here we represent the ratio = E{A^)/E{2^) for Sm and Gd iso- 
topes. We observe both in the theoretical and the experimental results that increasing the 
number of neutrons in the system the ratio varies from 2 -vibrational spectrum- 

to values close to ~ 3.3 that correspond to a rotational spectrum. This behavior is a 
signature of the presence of a shape transition in this region. We obtain a good agree- 
ment between the data and the theoretical results in the vibrational and rotational limits 
while some deviations are observed in the transitional isotopes N = 86 — 90. Precisely it 
is in these nuclei where the triaxial degree of freedom -neglected in this approximation- 
is more relevant. Furthermore, it is known that the octupole deformation -which violates 
parity symmetry- plays also an important role in these nuclei and could affect the qual- 
ity of the description in this region IfTTI . 

Finally, as it was mentioned above, the 150 Nd, 152 Sm and 154 Nd nuclei show experi- 
mentally excitation spectra similar to the one predicted by the X(5) solution E H [5J|. 
However, in the axial GCM calculations these isotopes possess a rotational character 
while the N = 88 nuclei are closer to the X(5) prediction -see horizontal dotted line in 
Fig. [5j This similarity is given not only in the ratio but also in the higher excited 
states and transition probabilities B(E2) Q (not shown). Nevertheless, these excitation 
spectra, though reproducing the X(5) spectrum, are extracted from intrinsic potentials 
that have not a clear connection to the X(5) potential (Fig. 151). Also Niksic and collab- 
orators fl9J obtain an excellent agreement with the experimental data for Nd with an 
intrinsic potential different to the one of the critical point. 

In summary, we have studied the spherical-prolate deformed shape transition in Nd, Sm 
and Gd isotopes with beyond mean field calculations using the Gogny D1S force. We 
have compared the results with the description of a first order shape-phase transition 
given by the collective hamiltonian showing that, contrary to this model, the triaxial de- 
gree of freedom plays an important role. Also the critical point potential is not observed. 
Both the vibrational and rotational empirical limits are obtained with GCM calculations 
although the agreement with the experimental data is worse in the transitional region, 
where some relevant degrees of freedom are neglected. Work is in progress to include 
them in future analyses. 
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